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\ ABSTRACT: We perform a detailed analysis of the non stationary solutions of the evolution 

(Fokker- Planck) equations associated to either stationary or non stationary quantum states by the 

4h 



stochastic mechanics. For the excited stationary states of quantum systems with singular velocity 
fields we explicitely discuss the exact solutions for the HO case. Moreover the possibility of modi- 
fying the original potentials in order to implement arbitrary evolutions ruled by these equations is 
discussed with respect to both possible models for quantum measurements and applications to the 
control of particle beams in accelerators. 



1. Introduction 

In a few papers [1] the analogy between diffusive classical systems and quantum systems has been 
reconsidered from the standpoint of the stochastic mechanics (SM) [2] , [3] , and particular attention 
was devoted there to the evolution of the classical systems associated to a quantum wave function 
when the conditions imposed by the stochastic variational principle are not satisfied (non extremal 
processes). The hypothesis that the evolving distribution converges in time toward the quantum 
distribution, constituted several years ago an important point in the answer by Bohm and Vigier 
to some criticisms to the assumptions of the Causal Interpretation of the Quantum Mechanics 
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(CIQM) [4]. In the quoted papers [1] it was pointed out that, while the right convergence was in 
fact achieved for a few quantum examples, these results could not be considered general as shown 
in some counterexamples: in fact not only for particular non stationary wave functions (as for a 
minimal uncertainty packet), but also for stationary states with nodes (namely with zeros) we do 
not seem to get the right asymptotic behaviour. For stationary states with nodes the problem is 
that the corresponding velocity field to consider in the Fokker-Planck equation shows singularities 
in the locations of the nodes of the wave function. These singularities effectively separate the 
available interval of the space variables into (probabilistically) non communicating sections which 
trap any amount of probability initially attributed and make the system non ergodic. 

In a more recent paper [5] it has been shown first of all that for transitive systems with 
stationary velocity fields (as, for example, a stationary state without nodes) we always have an 
exponential convergence to the right quantum probability distribution associated to the extremal 
process, even if we initially start from an arbitrary non extremal process. These results can also be 
extended to an arbitrary stationary state if we separately consider the process as confined in every 
configuration space region between two subsequent nodes. Moreover it has been remarked there 
that while the non extremal processes should be considered virtual, as trajectories in the classical 
Lagrangian mechanics, they can also be turned real if we modify the potential in a suitable way. 
The interest of this remark lies not only in the fact that non extremal processes are exactly what 
is lacking in quantum mechanics in order to interpret it as a totally classical stochastic process 
theory (for example in order to have a classical picture of a double slit experiment [6]), but also in 
the possibility of engineering some new controlled real evolutions of quantum states. In fact this 
could be useful to study (a) transitions between stationary states (b) possible models for measure 
theory [3] and (c) control of the particle beam dynamics in accelerators [7]. In a sense the SM is 
also a theory, independent from quantum mechanics, with applications in several physical fields, 
in particular for systems not perfectly described by the quantum formalism, but whose evolution 
is correctly controlled by quantum fluctuation: the so called mesoscopic or quantum-like systems. 
This behaviour characterizes, for example, the beam dynamics in particle accelerators and there 
is evidence that it could be described by the stochastic formalism of Nelson diffusions [1], [7]. 
Of course in this model trajectories and transition probabilities always are perfectly meaningful 
and, to study in detail the evolution of the probability distributions, and in particular to try to 
understand if and how it is possible to realize controlled evolutions, it is necessary to determine 
the fundamental solutions (transition probability densities) associated by SM to every quantum 
state in consideration: a problem dealt with in the following sections. 



2. Fokker-Planck equations for stochastic mechanics 

SM is a generalization of classical mechanics based on the theory of classical stochastic processes 
[2]. The variational principles of Lagrangian type provide a foundation for it, as for the classical 
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mechanics or the field theory [3]. In this scheme the deterministic trajectories of classical me- 
chanics are replaced by the random trajectories of diffusion processes in the configuration space. 
The surprisig feature is that programming equations derived from the stochastic version of the 
lagrangian principle are formally identical to the equations of a Madelung fluid [8], the hydrody- 
namical equivalent of the Schrodingcr equation in the Stochastic Interpretation of the Quantum 
Mechanics (SIQM) [9]. On this basis, it is possible to develop an interpretative scheme where the 
phenomenological predictions of SM coincide with that of quantum mechanics for all the experimen- 
tally measurable quantities. Within this interpretative code the SM is nothing but a quantization 
procedure, different from the ordinary ones only formally, but completely equivalent from the point 
of view of the physical consequences. Hence we consider here the SM as a probabilistic simulation 
of quantum mechanics, providing a bridge between this fundamental section of physics and the 
stochastic differential calculus. However it is well known that the most peculiar features of the in- 
volved stochastic processes, namely the transition probability densities, seem not always enter into 
this code scheme: in fact, if we want to check experimentally if the transition probabilities are the 
right ones for a given quantum state, we are obliged to perform repeated position measurements 
on the quantum system; but, according to quantum theory, the quantum state changes at every 
measurement (wave packet reduction), and since our transition probabilities are associated to a 
well defined wave function it will be in general practically impossible to experimentally observe 
a well defined transition probability. Several ways out of these difficulties have been explored: 
for example stochastic mechanic scheme could be modified by means of non constant diffusion 
coefficients [I]; or alternatively it would be possible to modify the stochastic evolution during 
the measurement [10]. Here we will rather assume that the processes which do not satisfy the 
stochastic variational principle still keep a physical meaning and that tey will rapidly converge (in 
time) toward the processes associated to quantum states. Indeed on the one hand any departure 
from the distributions of quantum mechanics will quickly be reabsorbed in the time evolution, at 
least in many meaningful cases; and on the other hand the non standard evolving distributions 
could be realized by suitable quantum systems for modified, time dependent potentials which may 
asymptotically in time rejoin the usual potentials. 

SM is a model intended to achieve a connection between quantum mechanics and classical 
random phenomena: here we will recall a few notions in order to fix the notation. The position of 
a classical particle is promoted to a vector Markov process £(i) defined on some probabilistic space 
(f2, P) and taking values in R 3 . We suppose that this process is characterized by a pdf /(r, i) 
and a transition pdf p(r, t\ r', t') and satisfies an ltd stochastic differential equation of the form 



where Vj are the components of the forward velocity field. However here Vj are not given a 
priori, but play the role of dynamical variables and are subsequently determined on the basis of a 
variational principle, namely on the basis of a dynamics. On the other hand rj(t) is a Brownian 




(2.1) 
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process independent of and such that 



F, t (dr)j(t)) =0, ^(drjji^drjkit)) =2DS jk dt 2.2 

where c?r/(t) = ry(t + dt) — ij(t) (for dt > 0), -D is a diffusion coefficient, and E t are the conditional 
expectations with respect to £(t). In what follows we will limit ourselves to the case of the one 
dimensional trajectories, so that the Markov processes £(t) considered will always take values in R. 
Moreover we will suppose for the time being that the forces acting on the particle will be defined 
by means of a time- independent potential V(x). A suitable definition of the Lagrangian and of the 
stochastic action functional for the system described by the dynamical variables / and v allows 
us to select, by means of the principle of stationarity of the action, the processes which reproduce 
the quantum mechanics [2], [3]. In fact, while the pdf f(x,t) of the process satisfies, as usual, the 
Forward Fokker-Planck (FP) equation associated to (2.1) 

d t f = Dd 2 J - 8 x (vf) = d x (Dd x f - vf) , (2.3) 

the following choice for the Lagrangian field 

£(x, t) = y v 2 (x, t) + mDd x v(x, t) - V(x) (2.4) 

enables us to define a stochastic action funcional 

ft! 

A= / E£(£(t),t) dt (2.5) 

J to 

which leads, through the stationarity condition SA — 0, to the equation 

dtS + + V _ 2mD 2 ^fl- = (2.6) 

2m V7 

involving a field S(x, t) defined as 

S(x, t) = - ^ E (C(£(s), s) | £(t) = x) ds + E (5i(^(ti)) | $(t) = x) (2.7) 

where S\{ ■ ) = S( - ,t\) is an arbitrary final condition. Now the relevant remark is that if R(x, t) = 
\/f(x, t), and if we define 

V-(x, t) = R(x, t) e lS(x ^/ R (2.8) 



the equation (2.6) takes the form 



and the complex wave function tp will satisfy the Schrodinger equation 



h 2 

ihd t i) = Hii) = - — dli> + Vii), (2.io) 



provided that the diffusionc coefficient be connected to the Planck constant by the relation 

D=A (2.11) 

2m v ; 

This trail leading from classical stochastic processes (plus a dynamics) to quantum mechanics 
can also be trod in the reverse way following the line of reasoning of the SIQM which, as it is well 
known, is formally ruled by the same differential equations as the SM. If we start from the (one 
dimensional) Schrodinger equation (2.10) with the Ansatz (2.8), and if we separate the real and 
the imaginary parts as usual in SIQM [8], the function f = R 2 = \ip\ 2 comes out to be a particular 
solution of a FP equation of the form (2.3) with constant diffusion coefficient (2.11) and forward 
velocity field 

v(x,t) = -d x S+^-d x (lnR 2 ). (2.12) 

On the other hand the explicit dependence of v on the form of R clearly indicates that to have a 
solution of (2.3) which makes quantum sense we must pick-up just one, suitable, particular solution. 
In fact the system is ruled not only by the FP equation (2.3), but also by the second, dynamical 
equation (2.9), the so-called Hamilton- Jacobi-Madelung (HJM) equation, deduced by separating 
the real and imaginary parts of (2.10) (see [8]). The analogy between (2.3) and a FP equation, 
which looks rather accidental in a purely SIQM context, is more than formal since, as we have 
briefly recalled, the SM shows how to recover both the equations (2.3) and (2.9) (and hence the 
Schrodinger equation (2.10)) in a purely classical, dynamical stochastic context. 



3. The eigenvalue problem for the FP equation 

Let us recall here (see for example [11]) a few generalities about the pdf's (probability density 
functions) f(x,t) solutions of a one-dimensional FP equation of the form 

d t f = d 2 x (Df) - d x (vf) = d x [d x (Df) - „/] (3.1) 

defined for x € [a, b] and t > to, when D(x) and v{x) are two time independent functions such that 
D(x) > 0, v(x) has no singularities in (a, b), and both are continuous and differentiablc functions. 
The conditions imposed on the probabilistic solutions are of course 

f(x,t) > 0, a<x<b,to<t, 

f b (3-2) 
/ f(x,t)dx = 1 , h <t, 

J a 

and from the form of (3.1) the second condition also takes the form 

[d x (Df)-vf] ab = 0, t <t. (3.3) 

Suitable initial conditions will be added to produce the required evolution: for example the tran- 
sition pdf p(x, t\x ,t ) will be selected by the initial condition 

lim f(x, t) = f(x, t+) = S(x - x ) . (3.4) 
t->i+ 



It is also possible to show by direct calculation that 

h(x) = N' 1 e~ f [D ' {x) ~ v ( x)]/D{x)dx , N = [ e ~f [D '( x) ~ v ( x)]/D{x)dx dx (3.5) 

J a 

is an invariant (time independent) solution of (3.1) satisfying the conditions (3.2). Remark however 
that (3.1) is not in the standard self-adjoint form [12]; but if we define the function g(x, t) by means 
of 

f{x,t) = ^/h{x)g{x,t) (3.6) 
it would be easy to show that g(x, t) obeys now an equation of the form 

dtg = Cg (3.7) 

where the operator L defined by 

dip(x) 



Lip = Tx 



p(x) 



dx 



q{x)<p(x) , (3.8) 



with 

p(x) = D(x) > 0, 



, . \D'(x) - v(x)} 2 \D'(x)-v(x)Y (3-9) 

q{x) = 4D(S) 2 ' 

is now self-adjoint. Then, by separating the variables by means of g(x, t) — r y{t)G(x) we have 
7(f) = c _At while G must be solution of a typical Sturm-Liouville problem associated to the 
equation 

CG{x) + XG(x) = (3.10) 

with the boundary conditions 

\D'(a) - v(a)] G(a) + 2D(a)G'(a) = , 

(3.11) 

[D'(b) - v(b)] G{b) + 2D{b)G\b) = . 
It easy to see that A = is always an eigenvalue for the problem (3.10) with (3.11), and that the 



corresponding eigenfunction is yjh(x) as defined from (3.5). 

For the differential problem (3.10) with (3.11) we have that [12] the simple eigenvalues A„ 
will constitute an infinite, increasing sequence and the corresponding eigenfunction G n {x) will 
have n simple zeros in (a, b). For us this means that A = 0, corresponding to the eigenfunction 
Gq(x) — y/h(x) which never vanishes in (a, 6), is the lowest eigenvalue and that all other eigenvalues 
are strictly positive. Moreover the cigenfunctions will constitute a complete orthonormal set of 
functions in L 2 ([a,6]) [13]. As a consequence the general solution of (3.1) with (3.2) will have the 
form 

oo 

f(x,t) = Y,^~ Xnt VW)Gn{x) (3.12) 

n=0 

with c = 1 for normalization (remember that A = 0). The coefficients c„ for a particular solution 
are selected by an initial condition 

f(x,t+) = f (x) (3.13) 
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and are then calculated from the orthonormality relations as 



In particular for the transition pdf we have from (3.4) that 

_ G n (x ) 

\AOo) 



b f Q ( x )^ldx. (3.14) 



(3.15) 



Since Ao = and A„ > for n > 1, the general solution (3.12) of (3.1) has a precise time evolution: 
all the exponential factors in (3.12) vanish with t — > +oo with the only exception of the term n = 
which is constant, so that exponentially fast we will always have 



lim f(x,t) = CQy // h(x)Go(x) — h(x) , (3.16) 
namely the general solution will always relax in time toward the invariant solution h(x). 



4. Stationary quantum states 

Let us consider now a Schrodinger equation (2.10) with a time-independent potential V(x) which 
gives rise to a purely discrete spectrum and bound, normalizable states, and let us use the following 
notations for stationary states, eigenvalues and eigenfunctions: 

<M*,t) =0„(x)e-* B "*/ fi 

h 2 ( 4J ) 

H<pn = 4>n + V< t>n = Kfa ■ 

Taking into account the relation (2.11) the previous eigenvalue equation can also be recast in the 
following form 

D<% = K . (4.2) 

For these stationary states the pdf is the time independent, real function 

f n {x) = \i> n {x,t)\ 2 = <t> 2 n {x), (4.3) 

and 

S{x,t) = -E n t, R(x,t) = (t> n {x), (4.4) 
so that for our state the velocity field is 

This means that now v n is time-independent and it presents singularities in the zeros (nodes) of the 
eigenfunction. Since the n-th eigenfunction of a quantum system with bound states has exactly n 
simple nodes [12] that we will indicate with x\,...,x n , the coefficients of the FP equation (2.3) are 
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not denned in these n points and we will be obliged to solve it in separate intervals by imposing the 
right boundary conditions connecting the different sections. In fact these singularities effectively 
separate the real axis in n + 1 sub-intervals with walls impenetrable to the probability current. 
Hence the process will not have an unique invariant measure and will never cross the boundaries 
fixed by the singularities of v(x): if we start in one of the intervals in which the axis is so divided 
we will always remain there [14]. 

As a consequence we must think the normalization integral (3.2) (with a = — oo and b = +oo) 
as the sum of n + 1 integrals over the sub-intervals [xk,Xk+i] with k = 0,1,..., n (where we 
understand, to unificate the notation, that x = — oo and x n+ i = +oo). Hence for n > 1 we will 
be obliged to solve the equation (2.3) in every interval [xk,Xk+i] by requiring that the integrals 

f fc+1 f(x,t)dx (4.6) 

be kept at a constant value for t > t : this value is not, in general, equal to one (only the sum of 
these n+1 integrals amounts to one) and, since the separate intervals can not communicate, it will 
be fixed by the choice of the initial conditions. Hence the boundary conditions associated to (2.3) 
require the conservation of the probability in [xk, Xk+i], namely the vanishing of the probability 
current at the end points of the interval: 

[ Dd *f- v f] Xk , Xk+1 =°> *^*o- ( 4 - 7 ) 

To have a particular solution we must moreover specify the initial conditions: in particular we will 
be interested in the transition pdf p(x, t\xQ, to), which is singled out by the initial condition (3.4), 
since [1] the asymptotic approximation in L 1 among solutions of (2.3) is ruled by the asymptotic 
behavior of p(x, t\xo,t ) through the Chapman-Kolmogorov equation 

/+oo 
p{x,t\y,t )f{y,t+)dy. (4.8) 
-oo 

It is clear at this point that in every interval [xk, Xk+i] (both finite or infinite) we can solve the 
equation (2.3) along the guidelines sketched in the section 3 by keeping in mind that in [xfc,Xfc+i] 
we already know the invariant, time-independent solution (pn( x ) ( or i mor e precisely, its restriction 
to the said interval) which is never zero in this interval with the exception of the extremes Xk and 
Xk+i- Hence, as we have seen in the general case, with the position 

f(x,t) = <j> n (x)g(x,t) (4.9) 

we can reduce (2.3) to the form 

d t g = C n g (4.10) 
where C n is now the self-adjoint operator defined on [xk, Xk+i] by 

-q n (x)ip(x) (4.11) 



CMx) = ^ 



p(x) 



d(p(x) 
dx 
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where we have now 

p(x)=D>0; qn{x ) = ^M + <M. (4.12) 

To solve (4.10) it is in general advisable to separate the variables, so that we immediately have 
j(t) = c~ At while G must be solution of the Sturm-Liouville problem associated to the equation 

C n G(x) + XG(x) = (4.13) 

with the boundary conditions 

[2DG'(x)-v n (x)G(x)] XkXk+i =0. (4.14) 

The general behaviour of the solutions obtained as expansions in the system of the eigenf unctions 
of (4.13) has already been discussed in section 3. In particular we deduce from (3.12) that for the 
stationary quantum states (more precisely, in every subinterval defined by two subsequent nodes) 
all the solutions of (2.3) always converge in time toward the right quantum solution \<j) n \ 2 : a general 
result not contained in the previous papers [1]. As a further consequence a quantum solution 
defined on the entire interval (— oo, +oo) will be stable under deviations from its initial condition. 



5. Harmonic oscillator 

To see in an explicit way how the pdf's of SM evolve, let us consider now in detail the particular 
example of a quantum harmonic oscillator (HO) characterized by the potential 

n*) = f"V. (5.1) 

It is well-known that its eigenvalues are 

E n = hcj ( n + J ) ; n = 0,l,2... (5.2) 



while, with the notation 



cr, 



2 



2 



h 



(5.3) 



2mw ' 
the eigenfuncions are 

K (x) = - l = e-* 2 / 4CT o Hn (^j=) (5.4) 



where H n are the Hermite polynomials. The corresponding velocity fields are easily calculated and 
are for example 

v {x) = —ux , 

v 1 {x) = 2— ujx, (5.6) 



v 2 (x) = Awal — j ~ uox ■ 

X <J Q 
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with singularities in the zeros Xk of the Hermite polynomials. If we now keep the form of the 
velocity fields fixed we can consider (2.3) as an ordinary FP equation for a diffusion process and 
solve it to see the approach to the equilibrium of the general solutions. When n = the equation 
(2.3) takes the form 

dtf = ua 2 dlf + uoxd x f + uf (5.7) 
and the fundamental solution comes out to be the Ornstein-Uhlenbeck transition pdf 

p(x,t\x ,t ) = -J—e-I*-* WW*) t {t > to) (5 . 8) 

where we used the notation 

a(t) = xoe'^-^ , a 2 (t) =a 2 [l-e^ t -*°)] , (t>t ). (5.9) 

The stationary Markov process associated to the transition pdf (5.8) is selected by the initial, 
invariant pdf 

f(x) = —^—c- x2 ' 2a « (5.10) 
(ToV 2tt 

which is also the asymptotic pdf for every other initial condition when the evolution is ruled by 
(5.7) (see [1]) so that the invariant distribution plays also the role of the limit distribution. Since 
this invariant pdf also coincides with the quantum stationary pdf </>q = IV'ol 2 the process associated 
by the SM to the ground state of a quantum HO is nothing but the stationary Ornstein-Uhlenbeck 
process. 

For n > 1 the solutions of (2.3) are no more so easy to find and, as discussed in the previus 
section, we will have to solve the eigenvalue problem (4.13) which, with e = h\, can be written as 

- £ G " {x) + (y^ 2 ~ ^ ^r~) G[x) = eG{x) ' (5 - n) 

in every interval [xfe,Xfc+i], with fc = 0,1, ... ,n, between two subsequent singularities of the v n 
field. The boundary conditions at the endpoints of these intervals, deduced from (4.7) through 
(4.9), are 

[4> n G' - 4>' n G\ Xk , Xk+1 = (5.12) 
and since <j) n (but not <j>' n ) vanishes in Xk,Xk+i, the conditions to impose are 

G(x k ) - G(x k+1 ) = (5.13) 
where it is understood that for xq and x n +i we respectively mean 

lim G(x)=0, lim G{x) = . (5.14) 



x — > — oo 



It is also useful at this point to give the eigenvalue problem in an adimensional form by using the 
new adimensional variable x/gq (which will still be called x) and the eigenvalue fi = \/u> = e/huj. 
In this way the equation (5.11) with the conditions (5.13) becomes 

... . fx 2 2n + l \ . . 

y"(x) u y(x) = 

W V 4 2 V yy ' (5.15) 

y{xk) = y{xk+i) = o 
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where x, Xk,Xk+i are now adimensional variables. If fi m and y m (x) are the eigenvalues and eigen- 
functions of (5.15), the general solution of the corresponding FP equation (2.3) will be 

f(x,t) = V c m e-^ t 4> n {x)y m (-) . (5.16) 

Of course the values of the coefficients c m will be fixed by the initial conditions and by the obvious 
requirements that f(x,t) must be non negative and normalized (on the whole x axis) along all its 
evolution. Two linearly independent solutions of (5.15) are 

V™ = e-^M (- £±!L t I , = xe -^ M (_ H+llzl, I , (5.17) 

where M(a,b;z) are the confluent hypergeometric functions. 

We consider first the case n = 1 (xo = — oo, x\ = and x 2 = +oo) so that (5.15) will have to 
be solved separately for x < and for x > with the boundary conditions y(0) = and 

lim y(x) = lim y(a;) = 0. (5.18) 

A long calculation [5] shows that the transition pdf is now 

x e -lx-a(t)] 2 /2* 2 (t) _ e -[x+a(t)] 2 /2a 2 (t) 

p(x, t\x , t ) = — — = (5.19) 

a(t) a(t)V2ir 

where a(t) and <r 2 (t) are defined in (5.9). It must be remarked however that (5.19) must be 
considered as restricted to x > when x > and to x < when x < 0, and that only on these 
intervals it is suitably normalized. In order to take into account at once both these possibilities 
we can also introduce the Heavyside function Q(x) so that for every x ^ we will have 

x e -lx- a (t)] 2 /2a 2 (t) _ c -[x+ a (t)] 2 /2* 2 (t) 

p{x,t\x ,t ) = e{xx ) — t — . (5.20) 

a{t) a(t)v2TT 

This completely solves the problem for n = 1 since from (4.8) we can now deduce also the evolution 
of every other initial pdf. In particular it can be shown that 

2 —x 2 /2a 2 

lim p(x,t\x ,t ) = 2O(xx ) X ^ " =29(xx )^(x), (5.21) 
and hence, if f(x, t$) = fo(x) is the initial pdf, we have for t > to 



/+oo 
p(x,t\y,t )fo(y)dy 
-oo 



/+oo 
&{xy)f {y)dy = T{q-x)<t>\{x), 
-OO 



(5.22) 



where we have defined the function 



T(q;x)=qQ(x) + (2-q)Q(-x); q = 2 f {y)dy. (5.23) 

Jo 
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Remark that when q = 1 (namely when the initial probability is equally shared on the two real 
semi-axis) we have F(l;x) = 1 and the asymptotical pdf coincides with the quantum stationary 
pdf <j>i (x); if on the other hand 5 ^ 1 the asymptotical pdf has the same shape of <j>\{x) but with 
different weights on the two semi-axis. 

If then n = 2 we have xq — —00, x\ = — 1, xi = 1 and X3 = +00, and the equation (5.15) 
must be solved in the three intervals (— 00,— 1], [—1,1] and [l,+oo), but the eigenvalues and 
eigenfunctions are now not easy to find so that a complete analysis of this case (and of every other 
case with n > 2) has still to be elaborated. At present only a few indications can be obtained 
numerically [5]: for example it can be shown that, beyond /i = 0, the first eigenvalues in the 
interval [—1,1] can be calculated as the first values such that 

M (-^±l,l')=0 (5.24) 



2 '2' 2, 

and are \i\ ~ 7.44, /i 2 ~ 37.06, /13 ~ 86.41. Also for the unbounded interval [l,+oo) (the analysis 
is similar for (—00, —1]) the eigenvalues are derivable only numerically. 



6. Controlled evolutions 

It is important to remark now that solutions of the type (5.8) and (5.20), and any other solution 
different from |</>„| 2 , are not associated to quantum mechanical states solutions of (2.10); in other 
words, they define processes that satisfy neither the stochastic variational principle [3] nor the 
Nelson dynamical equation [2]. That notwithstanding these processes still keep an interesting 
relation with the quantum mechanics. In fact to every solution f(x,t) of a FP equation (3.1), 
with a given v(x, t) and the constant diffusion coefficient (2.11), we can always associate the wave 
function of a quantum system if we take a suitable time-dependent potential. This means in 
practice that even the virtual (non optimal) processes discussed in this paper can be associated to 
proper quantum states, namely can be made optimal provided that the potential V(x) of (2.10) 
be modified in a new V(x,t) in order to control the evolution. 

Let us take a solution f(x,t) of the FP equation (3.1), with a given v(x,t) and a constant 
diffusion coefficient (3.3): if we define the functions R(x,t) and W(x, t) from 

f(x,t) = R 2 (x,t), v{x,t) = 8 x W(x,t), (6.1) 

if we remember from (2.12) that the following relation must hold 

mv = dxS + h^=d x S+^^=d x (s + ^ln/) (6.2) 

where / is an adimensional pdf (it is the argument of a logarithm) obtained by means of a suitable 
and arbitrary multiplicative constant, and if S(x, t) is supposed to be the phase of a wave function 
as in (2.8), we immediately get the equation 

S(x, t) = mW{x, t) - ^ In f(x, t) - 0{t) (6.3) 
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which allows us to determine S from / and v (namely W) up to an additive arbitrary function of 
the time 9 it). However, in order that the wave function (2.8) with the said R and S be a solution 
of a Schrodingcr equation, we must also be sure that the HJM equation (2.9) is satisfied. Since 
S and R are now fixed, the equation (2.9) must be considered as a relation defining a potential 
which, after a short calculation, becomes 

n 



Vix, t) = ^d 2 x In / + | (d t In / + vd x In /) - - md t W + 9. 



(6.4) 



Of course if we start with a quantum wave function for a given potential and if we pick up as a 
solution of (2.3) exactly / = R 2 , the formula (6.4) will correctly give back the initial potential, as 
can be seen for both the ground state and the first excited state of the HO which (by choosing 
respectively 9{t) = hut/2 and 0(f) = 3Uuit/2, which amounts to suitably fix the zero of the potential 
energy) give as result the usual harmonic potential (5.1). 

If on the other hand we consider for example the (non stationary) fundamental solution (5.8) 
associated to the velocity field Vo(x) of (5.6) for the case n = of the HO (we put to = to simplify 
the notation) we have already remarked that it does not correspond to a quantum wave function 
whatsoever. However a short calculation shows that, by choosing 

1 



0(t) 



~2~ 



2<7 2 



2 tanhejt 



and the time-dependent controlling potential 



Vix,t) = — 



ait) 



ait) 



9 9 

vrno x 



Tiuj 
~~2 



9 9 

muj x 



it -> +oo) , 



it -» +oo) 



(6.5) 



(6.6) 



we can define a quantum state (a wave function solution of a Schrodinger equation) which realizes 
the required evolution (5.8). Of course the fact that for t — ► +oo we recover the harmonic potential 
is associated to the fact, aleady remarked, that the usual quantum pdf (j^ix) is also the limit 
distribution for every initial condition and in particular also for the pdf (5.8). In the case n = 1, 
with viix) from (5.6) and the transition probability (5.20) as given non-stationary solution, the 
calculations are lenghtier. However if we define 



F(x,t) = 
xot{t) 



-[x- a (t)] 2 /2a 2 (t) 



ait)V2Tr 
xa(t) F{x,t) +G{x,t) 
a 2 it) F(x,t) - G(x,t) 



Gix,t) = 



[x+a(t)] 2 /2a 2 (t) 



er(t)V27f 

T{x) = 



tanh x 



and if we choose 

• hcj ( 4a 2 2a 2 a 2 (t) 
we have as time dependent potential for every x ^ 



it -> +oo) 



(6.7) 
(6.8) 

(6.9) 



Vix,t) 



2 2 



9 9 



24 



+ huj l - 

it -» +oo) . 



■ T 



xa\ 



4mx 2 



T B) 



(6.10) 
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In this case the asymptotic potential is the usual harmonic potential, but we must consider it 
separately on the positive and negative x semi-axis since in the point i = 0a singular behaviour 
would show up when t — * 0. This means that, also if asymptotically we recover the right potential, 
this will be associated with new boundary conditions in x = since we will be obliged to keep the 
system bounded on the positive (for example) semi-axis. 



7. Modelling transitions 

The explicit knowledge of the transition pdf of the type (5.8) and (5.20), And the possibility of 
turning optimal any suitable (/, v) state by a right choice of V(x, t) enable us also to explore the 
possibility of modelling evolutions leading, for example, from the pdf of a given stationary state to 
another (decays and excitations) . In fact a spontaneuous generalization of this idea hints to the 
possibility of modelling evolutions ffrom a given, arbitrary pdf and the pdf of an eigenfuncition of 
some observable: something which could become an element for very simple models of quantum 
measurements where we try to dynamically describe the wave packet collapse. As a first example 
let us consider the transition between the invariant pdf 's 

^ (7.1) 

f l{x ) = <Pl(x) = n ^=e-* 2 /^. 

If for instance we choose to describe the decay l^Owc should just use the Chapman-Kolmogorov 
equation (4.8) with (5.8) as transition pdf and fi(x) as initial pdf (to = 0). An elementary 
integration will show in this case that the resulting evolution takes the form 

./UoOM) = P 2 (t)Mx) +1 2 {t)h{x) (7.2) 

where we used the notation 

(3 2 (t) = 1 - c- 2wt , 7 (t) = e""* . (7.3) 

Taking now vq(x) from (5.6) and the evolving pdf from (7.2) and putting them in (6.4) (remark 
that, since w is stationary, d t W — 0) we get the following form of the controlling potential: 

9 9 
777 in x 

V(x,t) = — 2hwU(x/a ;p/ 7 ) (7.4) 

where 



In our example the parameter 



&2( t ) = ^ = e 2w * - 1 (7.6) 



7 2 (t) 
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is such that 6 2 (0 + ) = and 6 2 (+oo) = +00 and hence U goes everywhere to zero for t — ► +00, 
but is everywhere 1 with a negative singularity in x = for t — > + . As a consequence, while for 
t — > +00 the controlling potential (7.4) behaves like the HO potential (5.1), for £ > + it presents 
an unessential shift of —2tuv in the zero level, but shows also a deep negative singularity in x = 0. 

Apart from this singular behaviour of the controlling potential, a problem arises from the form 
of the phase funcion S. In fact from (6.3) we easily have for our decay 



S(x,t) = ~\n 



f3 2 {x,t) + ^ 2 {x,t) 



(7.7, 



so that in particular we have 

S(*,0+) = -§ln4, (7-8) 

while we would have expected that initially our phase function be independent from x as for every 
stationary wave function: this means that in our supposed evolution the phase function presents a 
discontinuous behaviour for t — > + . The problem arises here from the fact that in our simple model 
we initially have a stationary state characterized by a ddp /1 (a;) and a velocity field v\ (x) , and then 
suddenly, in order to start the decay, we suppose the same /1 embedded in a different velocity field 
vq{x) which drags it toward a new stationary fo{x). This discontinuous change from v\ to v is of 
course responsible for the remarked discontinuous change in the phase of the wave function. Hence 
a more realistic model for a controlled transition must take into account a continuous and smooth 
(albeit widely arbitrary) modification of the initial velocity field into the final one, a requirement 
which compels us to consider a new class of FP equations with time-dependent velocity field v(x,t). 
In particular to achieve the proposed controlled decay between two stationary states we should 
solve an evolution equation with a velocity field v(x,t) continuously, and possibly smoothingly, 
going from Vi(x) to vo(x); but this seems at present beyond the reach of our possibilities since 
every reasonable such v(x,t) field has proven intractable from the point of view of the solution of 
the FP equation (2.3). However we can show the results for another meaningful example which 
does not present the same technical difficulties of the decay between two stationary states: namely 
the controlled evolution from a coherent oscillating packet in a HO, and the ground state of the 
same HO. 

To do this we will recall a simple result [1] which indicates how to find the solutions of a 
particular class of evolution equations (2.3) which contains the situation of our proposed example. 
If the velocity field of the evolution equation (2.3) has the linear form 

v(x,t) = A(t) + B(t)x (7.9) 

with A(t) and B(t) continuous functions of time, then there are always solutions of the form 
Af((i(t), v{tj) where fi{t) and v(t) are calculated from the differential equations 

(j,'(t) - B{t)n{t) = A(t) ; v'{t) - 2B(t)u(t) = 2D (7.10) 
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with suitable initial conditions. On the other hand the (non stationary) wave function of the 
oscillating coherent wave packet with initial displacement a is 



i>c{x,t) = ( ) exp 



1/4 

J 



(x — a cos cut) 2 . (Aaxsinujt — a 2 sin2w< cut 



4a 2 ' \ 8a 2 



(7.11) 



so that the corresponding forward velocity field will be 

v c (x, t) = auo{cosujt — svciLut) — ujx , (7-12) 

namely it will have the required form (7.9) with A(t) = auj(cosujt — sinuit) and B(t) = — w, while 
the position pdf will be 

f c (x,t) - \^ c {x,t)\ 2 = f Q (x- a cos ut). (7.13) 

Now it is very easy to show that when B(t) = —u, as in the case of our wave packet, there are 
stable, coherent (non dispersive) solutions with v(t) = <Tq of the form N{n(t), ctq), namely of the 
form 

f(x,t) = f {x-»(t)). (7.14) 

Of course the time evolution of such coherent solutions can be determined in one step, without 
implementing the two steps procedure of first calculating the transition pdf and then, through the 
Chapman-Kolmogorov equation, the evolution of an arbitrary initial pdf. On the other hand if 
we compare (5.6) and (7.12) we see that the difference between vo and v c consists in the first, 
time dependent term of the second one; hence it is natural to consider the problem of solving the 
evolution equation (2.3) with a velocity field of the type 

v(x,t) = Ait) — ujx 

(7.15) 

A(t) = au;(cosu)t — smujt)F(t) 

where F(t) is an arbitrary function varying smoothly between 1 and 0, or vice verssa. In this case 
the evolution equation (2.3) still has stable, coherent (non dispersive) solutions of the form (7.14) 
with a fi(t) dependent on our choice of F(t) through (7.10). 

A completely smooth transition from the coherent, oscillating wave function (7.11) to the 
ground state </>o (5.4) of the HO can now be achieved for example by means of the following choice 
of the function F(t): 

N / AT\ 

F(t) - 1 - (1 - c- nt ) N = J2(-l) k+1 ( ,)e-" kt (7.16) 



fe=i 

where 



= ^^, Lu k = kQ; t>0, N>2. (7.17) 
r 

In fact this F(t) goes monotonically from F(0) = 1 to F(+oo) = with a flex point in r (which 
can be considered as the arbitrary instant of the transition) where its derivative F'(t) is negative 
and grows, in absolute value, logarithmically with N. The condition N > 2 also guarantees that 
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F'(0) = 0, and hence that the controlling potential V(x,t) of (6.4) will continuously start at t = 
from the HO potential (5.1), and eventually come back to it for t — > +oo. Finally the phase 
function S(x, t) too will change continuously from that of tp c to that of the HO ground state. A 
long but simple calculation will now show that the explicit form of the controlling potential is 

r 2 N / N\ 

V{x,t)=muj 2 —-mLjaxy2(-l) k+1 ( ,)[U k {t)LO k e-" kt - Wkuoe^] (7.18) 

2 fc =i W 

where 

. . . 2lo 2 sin tot — wi cos cot 
Uk(t) = smut H 



(Cjfe - W) 2 + LO 2 



w k = i + ( 2 " 2 7j 2 =^ fc (f 



(7.19) 



The parameters r and iV, with the limitations (7.17), are free and connected to the particular form 
of the transition that we want to implement. We conclude this section by remarking that, in a 
HO, the transition between a coherent, oscillating wave packet and the ground state is a transition 
between a (Poisson) superposition of all the energy eigenstates to just one energy eigenstate: an 
outcome which is similar to that of an energy measurement, but for the important fact that here 
the result (the energy eigenstate) is deterministically controlled by a time dependent potential. In 
fact our controlled transition does not produce mixtures, but pure states (eigenstates) and in some 
way realizes a dynamical model for one of the branches of a measurement leading to an eigenvalue 
and an eigenstate. 



8. Beam dynamics in particle accelerators 

As a model which tries to put in evidence the classical aspects of the quantum physics, the SM 
seems especially suitable to the description of systems whose nature in some sense lies between 
classical and quantum: the so called mesoscopic or quantum- like systems [15]. We will propose now 
a few preliminary remarks about the possibility of making use of this characteristic in a particular 
physical domain [7] . The dynamical evolution of beams in particles accelerators is a typical example 
of mesoscopic behaviour. Since they are governed by external electromagnetic forces and by the 
interaction of the beam particles among themselves and with the environment, charged beams 
are higly nonlinear dynamical systems, and most of the studies on colliding beams rely either on 
classical phenomena such as nonlinear resonances, or on isolated sources of unstable behaviors 
as building blocks of more complicated chaotic instabilities. This line of inquiry has produced 
a general qualitative picture of dynamical processes in particle accelerators at the classical level. 
However, the coherent oscillations of the beam density and profile require, to be explained, some 
mechanism of local correlation and loss of statistical independence. This fundamental observation 
points towards the need to take into account all the interactions as a whole. Moreover, the overall 
interactions between charged particles and machine elements are really nonclassical in the sense 
that of the many sources of noise that are present, almost all are mediated by fundamental quantum 
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processes of emission and absorbtion of photons. Therefore the equations describing these processes 
must be, in principle, quantum. 

Starting from the above considerations, two different approaches to the classical collective 
dynamics of charged beams have been developed, one relying on the FP equation [16] for the beam 
density, another based on a mathematical coarse graining of Vlasov equation leading to a quantum- 
like Schrodingcr equation, with a thermal unit of emittance playing the role of Planck constant [17]. 
The study of statistical effects on the dynamics of electron (positron) colliding beams by the FP 
equation has led to several interesting results, and has become an established reference in treating 
the sources of noise and dissipation in particle accelerators by standard classical probabilistic 
techniques [18]. Concerning the relevance of the quantum-like approach, at this stage we only 
want to point out that some recent experiments on confined classical systems subject to particular 
phase-space boundary conditions seem to to be well explained by a quantum-like (Schrodingcr 
equation) formalism [19]. In this approach [20] the (one dimensional) transverse density profile of 
the beam is described in terms of a complex function, called beam wave function, whose squared 
modulus give the transverse density profile of the beam. This beam wave function satisfies a 
Schrodinger-like equation where h is replaced by the transverse beam emittance e: 



On the other hand a recently proposed model for the description of collective beam dynamics in 
the semiclassical regime [21] relies on the idea of simulating semiclassical corrections to classical 
dynamics by suitable classical stochastic fluctuations with long range coherent correlations, whose 
scale is ruled by Planck constant. This elaborates a hypothesis first proposed by Calogero [22] in his 
attempt to prove that quantum mechanics might be interpreted as a tiny chaotic component of the 
individual particles' motion in a gravitationally interacting universe. The virtue of the proposed 
semiclassical model is twofold: on the one hand it can be formulated both in a probabilistic FP 
fashion and in a quantum-like (Schrodinger) setting, thus bridging the formal gap between the two 
approaches. On the other hand it goes further by describing collective effects beyond the classical 
regime due to the semiclassical quantum corrections. 

Since we are interested in the description of the stability regime, when thermal dissipative 
effects are balanced on average by the RF energy pumping, and the overall dynamics is conser- 
vative and time-reversal invariant in the mean, the choice to model the random kinematics with 
the Nelson diffusions, that are nondissipative and time-reversal invariant, is particularly natural. 
The diffusion process describes the effective motion at the mesoscopic level (interplay of thermal 
equilibrium, classical mechanical stability, and fundamental quantum noise) and therefore the dif- 
fusion coefficient is set to be the semiclassical unit of emittance provided by qualitative dimensional 
analysis. In other words, we simulate the quantum corrections to classical deterministic motion 
(at leading order in Planck constant) with a suitably defined random kinematics replacing the 
classical deterministic trajectories. Therefore, apart from the different objects involved (beam 
spatial density versus Born probability density; Planck constant versus emittance), the dynamical 



dip(x, z) 
dz 



e 2 d 2 %l>(x, z) 
~2 dx~ 2 



+ U(x, z)i[)(x, z) . 



(8.1) 
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equations of our model formally reproduce the equations of the Madelung fluid (hydrodynamic) 
representation of quantum mechanics. In this sense, the present scheme allows for a quantum-like 
formulation equivalent to the probabilistic one. 

With a few changes in the notation we can now reproduce, for the beam dynamics, the SM 
approach sketched in section 2. Let q(t) be the process representing some collective degree of 
freedom of the beam with a pdf p{x,t). Then, in suitable units, the basic stochastic kinematical 
relation is an Ito stochastic differential equation of the type (2.1) where the emittance e of the 
beam plays the role of a diffusion coefficient. Since we are interested in the stability regime of 
the bunch oscillations, the bunch itself can be considered in a quasi-stationary state, during which 
the energy lost by dissipation is regained in the RF cavities. In such a quasi-stationary regime 
the dynamics is, on average, invariant for time-reversal and we can define a classical effective 
Lagrangian L(q, q) of the system, where the classical deterministic kinematics is replaced by the 
random diffusive kinematics (2.1). The equations for the dynamics can then be obtained from the 
classical Lagrangian by means of the stochastic variational principles. 

Introducing now the time-like coordinate s = ct we get now the analog of the equations (2.3) 
and (2.6) in the form of a HJM equation 

d s S+^-2e 2 ^ + V(x,s) = 0, (8.2) 

and of a continuity equation 

d s p=-d x (pv). (8.3) 

Remark that now the symbol v no more represents the forward velocity fields, but rather the drift 
velocity connected to the forward and backward velocities by the relation 2v = +vr_\, and to 
the phase function by the relation v — d x S. The observable structure is now quite clear: Ew is the 
average velocity of the bunch center oscillating along the transverse direction; Eq gives the average 
coordinate of the bunch center; finally the second moment (Aq) 2 = E(q — E(q)) 2 determines the 
dispersion (spreading) of the bunch. The coupled equations of dynamics may now be used to 
achieve a controlled coherence: given a desired state (p, v) the equations of motion (8.2) and (8.3) 
can be solved to calculate the external controlling potential V(x, s) that realizes this state. 

General techniques to obtain localized quantum wavepackets as dynamically controlled systems 
in SM have already been introduced [23] . In this way one can construct for general systems either 
coherent packets following the classical trajectories with constant dispersion, or coherent packets 
following the classical trajectories with time-dependent, but at any time bounded dispersion. These 
results can now be extended also to the quantum-like description of the transverse dynamics of a 
particle beam and hence it will be possible to select a current velocity, by fixing the characteristics 
of the motion of the packet center, to determine the corresponding solutions of the FP (continuity) 
equation and finally to use the HJM equation as a constraint giving us the controlling device. The 
formal details of this program will be developed in a subsequent paper. 
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9. Concluding remarks 



It has been observed that the inverse problem of determining a controlling potential for a given 
quantum evolution in fact does not need to be formulated in terms of SM. Given two quantum 
wave function ipi and tp2 we could indeed design a new wave function vp(x,t), evolving from ipi to 
ip2 phig it, as required evolution, directly in the Schrodinger equation (2.10) and eventually deduce 
from that the form of the controlling potential. At first glance this seems to completely circumwent 
the need for a model like the SM: given an arbitrary evolving state we can always calculate the 
potential producing it. However about that two remarks are in order. 

First of all, from a purely technical point of view, the simplification introduced by this proce- 
dure shows up to be elusive. In fact we must remember that a quantum wave function has complex 
values and hence, if we simply take an arbitrary evolution, the resulting potential calculated from 
the Schrodinger equation (2.10) will also be complex. This means that, to have a real valued 
potential, we must impose some conditions on the supposed evolution. These conditions of course 
depend on the hypothesized form of ip. For example, if we fix the evolution of its modulus, the 
said condition will materialize in a partial differential equation on the phase function S of the wave 
function. On the other hand the use of the HJM equation (2.9) as the tool to solve the inverse 
problem always give a real valued potential as a result. 

However both the two proposed procedures are possible and, to identically posed questions, 
they will give identical answers. Given this obvius equivalence, the second remark is that our 
choice of the procedure will be operated on the basis of opportunity considerations. In both cases 
the result will be influenced by the starting hypothesis on the supposed evolution of the state ip 
modelling the transition from ipi to V>2- But, since the observable part of the wave function is its 
square modulus, namely the position pdf, the relevant hypothesis will be on its evolution. The 
phase function, or, equivalently, the velocity fields, are not directly observable, and hence are at 
first sight of secondary concern. Their importance become apparent only when we require that 
the potential be real or that the transitions show a realistic, smooth behaviour. Hence, depending 
on the specific problem we are dealing with, it could be more suitable to approach it in terms of a 
state given through a wave function tp, or in terms of a state given through the couple (/, v). The 
two approaches are certainly equivalent, but one may prove to be more suggestive. In particular 
that based on the SM equations seems to be better for the treatment of systems, like as the 
mesoscopic, quantum-like ones, which are well described by classical probabilistic models in terms 
of real space-time trajectories. 
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